Scalar phantom energy as a cosmological dynamical system 



O 
O 

(N 

a 

i> 

(N 

> 
o 
in 
m 
o 
o 
in 
o 

^ 

a: 

6 

03 



13 



L. Arturo Urena-Lope^ 

Institute) de Fisica de la Universidad de Guanajuato, 
A. P. E-143, CP. 37150, Leon, Guanajuato, Mexico 
(Dated: February 5, 2008) 

Phantom energy can be visualized as a scalar field with a (non-canonical) negative kinetic energy 
term. We use the dynamical system formalism to study the attractor behavior of a cosmological 
model containing a phantom scalar field <j> endowed with an exponential potential of the form 
V(<j>) = Vo exp(—\K<f>), and a perfect fluid with constant equation of state 7; the latter can be of the 
phantom type too. As in the canonical case, three characteristic solutions can be identified. The 
scaling solution exists but is either unstable or of no physical interest. Thus, there are only two 
stable critical points which are of physical interest, corresponding to the perfect fluid and scalar 
field dominated solutions, respectively. The most interesting case arises for > 7 > -A 2 /3, which 
allows the coexistence of the three solutions. The main features of each solution are discussed in 
turn. 

PACS numbers: 98.80.-k,95.35.+d,05.45.-a 



I. INTRODUCTION 

For many years, scalar fields have played the role of 
work-horses for many phenomenological models in mod- 
ern Cosmology. And the main issue is to find which so- 
lutions of the field equations are physically important to 
explain the universe we live in. 

The usual route was to get expertise in solving the 
coupled Einstein-Klein-Gordon (EKG) set of nonlinear 
differential equations, as it can be seen in the early liter- 
ature on scalar fields. However, it was recently realized 
that for a exponential potential the EKG equations could 
be seen as a plane autonomous system of equations 0. 
That is, the EKG equations could be seen as a dynami- 
cal system. This case became one more among the many 
examples of dynamical systems in Cosmology, see@- 

Even though the exponential potential is one of the 
most studied cases in the cosmology of scalar fields 
(see0, 0, |M and references therein), the way it is 
handled injlj allows a better understanding of the evo- 
lution of a universe with exponential potentials. The 
idea has been exploited in some papers that have used 
it successfully for diverse situations. For instance, in 
the cases of negative exponentials^, an exponential po- 
tential in general Robertson- Walker spacetimes[^, and 
multi-exponential potentials'^!; see also[Tl|- 

On the other hand, the idea of a scalar field with 
non-canonical kinetic energy has drawn the attention 
of cosmologists because such a scalar field can have, 
in princi ple , what is called a phantom equation of 
state[H [H Q E E H3, for which p/p < -1, be- 
ing p and p the pressure and the energy density of the 
phantom fluid, respectively. This new type of fluids vio- 
lates the so-called dominant energy condition, which for 
a perfect fluid is written as \p\ < p\ for a thoroughly 
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discussion on this condition and the stability of vacuum 
see 0,0- Thus, the case of phantom scalar fields should 
be taken with care, as some even think that a cosmolog- 
ical constant may be easier to justify 0, 0, The 
presence of a kind of phantom matter may be supported 
by supernovae type la observations, as phantom mat- 
ter seems to fit them better than a cosmological con- 
stant or a quintessence field 0, 0, |2jj. Other alter- 
natives, though, are a canonical scalar field climbing 
up its scal ar p otentialpH 122^ . quantum corrections on 
large scales [23l |24| . and a non-canonical complex scalar 
field |25|: mechanisms that may result in a phantom-like 
equation of state. 

Recently, the attractor behavior of phantom scalar 
field models, which have a negative kinetic term, was 
studied in 26] using the Hamilton-Jacobi formalism. 
Such study was restricted to cosmological models in 
which the phantom scalar field was the only matter 
present, using three different types of scalar potentials 
(power-law, exponential and cosine). 

In the present paper, our aim is to perform a more 
complete study of the attractor properties of a phan- 
tom scalar field endowed with a positive exponential po- 
tential, using the dynamical system formalism developed 
in[lj for canonical scalar fields. This study will include a 
perfect fluid with a constant barotropic equation of state. 

In this respect, it should be said that a similar anal- 
ysis was performed in|27l l28j . where a coupling between 
the phantom scalar field and a dust fluid was included. 
The main difference with that work is that we will not 
consider any coupling, and will not restrict the perfect 
fluid equation of state to the usual values; we will allow 
the possibility of super-stiff and phantom perfect fluids. 

There are many other authors that have studied the 
cosmologies containing a phantom field, and have ob- 
tained explicit solutions together with careful stability 
analysis, see [H E |H E [H HJ HI H3 and the 
long list of references there in. We will comment on any 
similarity or difference that may appear with respect to 
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those papers. 

A brief summary of the contents of this paper is as fol- 
lows. In Seem we will review the mathematics needed 
for a cosmological model with a perfect fluid and a phan- 
tom scalar field. In Sec. IIIII the critical points of the 
system and their stability are presented. The discussion 
is focused in the three different cases that appear accord- 
ing to the values the perfect-fluid equation of state can 
take. Finally, Sec. IIVI is devoted to conclusions. 



II. MATHEMATICAL BACKGROUND 

For simplicity, we shall restrict ourselves to the case 
of a homogeneous and isotropic universe with zero cur- 
vature, its metric given by the flat Robertson- Walker 
one|39j. In this universe, there is a perfect fluid with 
pressure p 1 and energy density p 1 related through a con- 
stant equation of state 7 in the form p 1 = (7 — l)p 1 . 
There is also a phantom scalar field <f> minimally cou- 
pled to gravity, whose Lagrangian density reads C = 
(l/2)9 M 09 AI (/>— V{4>), where V(<fr) is called the scalar po- 
tential. Therefore, the equations of motion are the known 
EKG equations 



H = —j [IPy- 



• dV 

-3H(/)+ — . 



(la) 
(lb) 

p 7 = -3#7p 7 , (lc) 
together with the constraint (Friedmann) equation 

1 • 



< iJ H;) 2 4M* 2 1 



(2) 



where dots mean derivative with respect to cosmic time, 
k 2 = 8ttG, H is the Hubble parameter, and a(t) is the 
scale factor of the universe. The above equations are the 
same as obtained in the canonical case except for the 
minus sign in front of the scalar kinetic term if) 2 . 

We will study the case of a scalar exponential potential 
of the form V(<f>) — Vo cxp(~ Xncj)) , where Vo and A are 
free parameters. The similarity with the canonical case 
suggests that we should make use of the dimensionless 
variables pj 



K(f) 

7^' 



y 



kVV 
V3H 



(3) 



and then Eqs. can be written as a plane-autonomous 
system 



-3a; + - [7 (l + x 2 - y 2 



2x 2 }x-^\y\i & ) 
[7 {l + x 2 -y 2 )- 2x 2 ]y-^Xxy, (4b) 



where now a prime denotes derivative with respect to the 
logarithm of the scale factor N = In a. The constraint 
equation JSJ) now reads 



1 = 



K 2 p<y 

3H 2 



(5) 



There are two scalar quantities we shall be interested 
in, which are the phantom density parameter and the 
effective phantom equation of state 7^, given by 



ilrf, = —x 2 



V 



= ~2x z 



(6) 



We notice that the scalar equation of state 7^ is always 
negative. 

It is interesting to note that there will be three different 
regions in the phase plane, characterized by 0^ < 0, < 
Q<t> < 1, and £1$ > 1, respectively. The boundary lines 
between the above regions are given by the hyperbola 
y 2 — x 2 = 1 and the straight lines y — ±x, and are shown 
with dashed curves in the figures below. 

As we will focus our attention in p 7 > 0, we see the 
Friedmann constraint @ requires all realistic trajectories 
in the phase plane {x,y} to comply with < 1. How- 
ever, in opposition to the canonical case, there is no rea- 
son to constraint only to positive values. As we shall 
see later, the existence of perfect-fluid supra- dominated 
eras cannot be discarded. 



III. CRITICAL POINTS AND STABILITY 

Eqs. (0J are written as an autonomous phase system 
of the form x' = f(x), and its so-called critical points 
Xo are solutions of the system of equations f(xo) = 0. 
To determine their stability we need to perform linear 
perturbations around the critical points in the form x = 
xo + u, which results in the following equations of motion 
u' = A4u, where 



dxA 



(7) 



A non-trivial critical point is called stable (unsta- 
ble) whenever the eigenvalues rrii of M. are such that 
Re(mi) < (Re(mi) > 0). If neither of the aforemen- 
tioned cases is accomplished, the critical point is called 
a saddle point Q- 

We applied the above procedure to Eqs. I@J. The 
critical points together with the stability analysis, and 
the corresponding values of Q<p and 7^ are shown in Ta- 
ble U 40] . For completeness, the eigenvalues of matrix Q) 
for the critical points A, B, C and D are shown in Ta- 
ble [Q] For solution E, the eigenvalues are 



m± 



7l 



(2- 7 ) ± / 87 (A 2 + 37) 



|2-7l 



A 2 (2 - 7) 



(8) 
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TABLE I: Critical points for a phantom scalar field endowed with an exponential potential. 



Label 


Pi 


xo 




Existence 


Stability 




70 


A 


PA 


(0,0) 




VA 


Saddle 





Undefined 


B 


PA 




xye) 


VA 


Stable 


1 


-A 2 /3 


C 


> 


(0,0) 




VA, V7 


Unstable for 7 > 2 
Saddle for 2 > 7 > 
Stable for 7 < 





Undefined 


D 


> 


(-A/V6 J% /1 + 


Ave) 


VA, V7 


Stable for 7 > -A 2 /3 
Saddle for 7 < -A 2 /3 


1 


-A 2 /3 


E 


> 


(73727/^37(7 


-2)/2A») 


VA, 7 > 2 
VA, 7 < 


Saddle for 7 > 2 
Saddle for > 7 > -A 2 /3 
Stable for 7 < -A 2 /3 


< 

-37/A 2 

> 1 
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TABLE II: Eigenvalues found for the critical points A, B, C 
and D shown in Table Q see Eq. J7J. 



Label 


m_ 


m+ 


A 


-3 





B 


-3- A 2 /2 


-A 2 


C 
D 


-3 + 37/2 
-3- A 2 /2 


3 7 /2 
-37 - A 2 



There are many similarities in the critical points with 
respect to the canonical case, so we briefly review the 
main solutions in the latter, see Table I in[l|. There are 
three solutions that are relevant for Cosmology, which are 
respectively dubbed the perfect fluid dominated solution, 
the scalar field dominated solution, and the scaling so- 
lution (for an interesting discussion on the existence of 
scaling solutions, seejsH). 

The fluid dominated solution is in general unstable, 
except in the case of a cosmological constant. The scalar 
field dominated solution only exists for A 2 < 6, in which 
case is the only late-time attractor and is inflationary 
for A 2 < 2. Also, this solution is a saddle point when a 
cosmological constant is present. 

For both cases 37 < A 2 and 6 < A 2 , the only stable so- 
lution is the scaling one, in which the scalar energy den- 
sity remains proportional to that of the perfect fluid with 
^0 = 37/A 2 . This last property has been widely used in 
some models of scalar field dark matter [3a. l37 | |41 | . 

As we shall see, the minus sign of the kinetic energy 
introduces relevant changes for the phantom case. There 
are three different situations we are about to revise now. 



A. 7 = 

As it was found in|27ll28l|. in a universe with a cosmo- 
logical constant p 7 = p\ and a phantom scalar field there 
is only one stable critical point, which is the scalar field 
dominated solution B. This the same as in the canoni- 
cal case, except that now there is not restriction on the 
values that A can take [53, [2^| . As we mentioned before, 



this also implies that an exponential potential will come 
to dominate over a cosmological constant too; this is to 
be expected since 7^ < 0. 

1. Only phantom scalar field 

It should be realized that the case in which no other 
matter except the phantom scalar field is present is not 
equivalent to the above case with 7 = 0. In the former, 
Eqs. (@J, and © can be fused into the single equation 

x' = -(zx+^-\^{l + x 2 ) . (9) 

Clearly, the only physically reasonable critical point is 
the same as point B 29j|42], which is again stable. This 
time, however, point A cannot exist. 

B. 7 > 

On the other hand, when a (7 ^ 0) perfect fluid is 
present, things become more interesting. For 7 > 0, the 
scaling solution E either does not exist or is not of phys- 
ical interest, and the perfect fluid dominated solution 
C is in general unstable. Therefore, we conclude that 
the (phantom) scalar field dominated solutions B and 
D are unavoidable if 7 > 0, which is the case for ordi- 
nary matter (radiation, dust, and even quintessence fields 
and a cosmological constant) 28] . The effective equation 
of state 70 is of the phantom type, and the scalar field 
climbs up its scalar potential. This is shown in Fig^for 

7=1; 

This is the first example in this work where we see 
supra- dominance of the perfect fluid over the phantom 
field, i.e. f2 7 > 1. Such a phenomenon, that passed 
unnoticed in previous published works |28ll29|. also arises 
if the perfect fluid and the phantom field are coupled, as 
in[13. 

Also in Fig. ^ we plotted trajectories that start with 
Q,<j, — (at the boundary lines y = ±x). We believe that 
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FIG. 1: The phase plane for 7=1 and A = 3. For 7 > 0, the 
only late-time attractor is the stable critical point D. Notice 
that, in opposition to the canonical scalar field, there is not re- 
striction on the value of A. We also show the phantom-divide 
(dashed) circumference, and some trajectories where the per- 
fect fluid is supra-dominant, see text for details. For com- 
pleteness, we show a representative trajectory with Q^, > 1 
(dotted curves), which never enters the physically interesting 
region < 1. For the case in which there is not a perfect 
fluid present, all trajectories lie on the hyperbola y 2 — x 2 — 1 
and finish at the critical point D. 
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FIG. 2: The phase plane for 7 = —4 and A = 3. For this case, 
the only late-time attractor is the stable critical point C, i.e., 
the (phantom) perfect fluid dominates over the scalar field, 
even if a stage of supra-dominance is allowed initially. For 
completeness, we show numerical solutions of the case Q$ > 
1 (dotted curves), for which point E is an attractor stable 
solution, which is in agreement with the stability analysis in 
Eq. ©. 



C. 7 < 



all of such trajectories are indeed part of longer trajecto- 
ries which may have started in the negative-O^, regions. 

It is also worth noticing that trajectories can start 
with f2^ < 0, but always end up with < fi^ < 1, 
but not vice-versa; this is because there are not stable 
critical points with 17^ < (in the case of coupling be- 
tween matter and the phantom field, there are stable 
critical points in which the perfect fluid matter is supra- 
dominant, see[27|]). Moreover, the parabola y 2 — x 2 = 1 is 
not traversable at all, as trajectories with Q<p > 1 never 
enters into the region < < 1 in order to reach the 
stable critical point C. 

To better visualize the evolution of the universe along 
the trajectories in Fig. ^ it is convenient to define an 
effective equation of state of the universe 7 ff, which for 
our case reads 

7eff = ■ =7 1 - y - x . (10) 

+ L V 7 J . 

The value j c s = defines, in general, an ellipse (hyper- 
bola) in the xy plane if0<7<2(7>2). That is to say, 
this ellipse (and this hyperbola) marks the phantom di- 
vide of the universe in the phase plane: j c s > (7 e g < 0) 
inside (outside) the ellipse (and the hyperbola). 

For 7=1, the phantom divide is a unitary circumfer- 
ence, as shown in Fig. Q A curious point is that some 
trajectories cross over the phantom divide twice. Even 
though the perfect fluid can be supra-dominant, the ef- 
fective equation of the universe rarely coincides with that 
of the perfect fluid. 



There is an important change in the analysis above if 
we allow for 7 < 0, that is, if the perfect fluid is also 
phantom. For 7 < — A 2 /3, D is a saddle point (the scalar 
field dominated solution cannot be sustained) , and point 
E is again of no physical interest even if it is a stable node. 
The only possibility is then the (phantom) perfect fluid 
dominated solution C (the perfect fluid is more phantom 
than the scalar field) . The phase plane for an instance of 
this case is shown in Fig. [21 

There is still one extra case we have to care about. 
For > 7 > — A 2 /3, both points C and D exist and are 
stable. Moreover, point E is also physically allowed, but 
this time is a saddle point. The numerical experience, as 
shown in Fig. [31 indicates that some of the trajectories 
pass nearby point E before finishing at one of the late- 
attractors C and D. How close a trajectory is to point E 
depends on the initial conditions Xo. 

We also see that all trajectories with initial perfect 
fluid supra-dominance never reach point D, and then al- 
ways end up at critical point C. As for the O > 1 case, 
we see that all trajectories end up at point D, but again 
without crossing over the < 1 region. 

As a final note, it can be seen that the phantom-divide 
curve for 7 < is a hyperbola, see Eq. I|10fl ; and then "f e g 
is positive only well within the (O^ > l)-region of the 
phase plane. 

IV. CONCLUSIONS 

The analysis of the EKG equations from the point of 
view of dynamical systems is very fruitful, as the appear- 
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FIG. 3: The phase plane for 7 = — 1 and A = 3. For this case, 
there are two late-time attractors, which are the critical points 
C and D; though the saddle point E is also present. Which 
of the stable critical points is reached at the end, depends on 
the initial values xo. Notice that point D is also a late-time 
attractor for Q$ > 1 (dotted curves). 

ance of critical points helps a lot to reveal the role scalar 
fields play in many cosmological models. 

We have studied the case of a phantom scalar field 
endowed with an exponential potential, and found that 
there are in general only two stable solutions which are of 
physical interest. For matter with 7 > 0, it is confirmed 
that the scalar field dominated solution is the only late- 
time attractor. 

There is for this case the same problem as in the canon- 
ical scalar field 0. In order to have a cosmological evolu- 
tion with standard radiation and matter dominated eras, 
the scalar field contribution should be almost negligible 
at early times. If the phantom scalar field were initially 
present at the time of inflation, the inflaton field could 
not avoid a phantom era unless f2^ C 1 at the begin- 
ning of inflation. This is because the critical point A is 
a saddle point, as we have numerically confirmed, but 
trajectories move away from it only linearly with N. 

As in the canonical case, once the radiation and matter 
dominated eras appear, the phantom scalar field density 
parameter grows like ~ a 37 , which also means that 
the scalar energy density remains almost constant until it 
reaches its late-time attractor solution. Hence, in order 
to avoid a disturbing phantom dominated universe, we 
need the same constraint that appears in the standard 
case, that is, the initial scalar density parameter should 
be such that < Pjfi/pjj, i n order not to become the 
dominant component at time to. 



The novel case that appears for a phantom equation 
of state 7 is the coexistence of three critical points if > 
7 > — A 2 /3; a case that can also happen when matters are 
coupled 27] . Allowing for more than one phantom type 
of matter leads to a non-trivial evolution; we learn in 
this case that the ultimate phantom regime will depend 
on the kind of phantoms present. For instance, it is seen 
from Fig.[3]that the most phantom matter may not be the 
dominant one at the end. This result may be relevant in 
models in which more than one type of phantom matter 
coexist together. In[34|], the authors report a lower bound 
in the effective equation of state in the case of many 
phantom scalars with an O(N) symmetry. This suggests 
that the classical solutions of phantom scalar fields may 
have hidden properties worth investigating. 

Another result is that supra- dominance of the perfect 
fluid 7 > I cannot be discarded in the case of a phan- 
tom scalar field, because the scalar field contribution fi^ 
cannot be, in principle, bounded from below. If phantom 
scalar fields were permitted in Nature, there would have 
to be a mechanism that can prevent to become too 
negative. 

Strangely enough, for < 7 < 2 the effective equation 
of state of the universe jeg can be negative even in the 
case of supra-dominance; moreover, it can change from 
negative to positive values, and vice-versa, along some 
trajectories. 

The model presented in this work has the same prop- 
erties that appears in any single field-phantom model in 
which 70 < 0: big rip singularity and no crossing of the 
70 = barrier, see[3(J for more details. 

In studying such properties, the dynamical system for- 
malism is not only useful but also easy. Apart from help- 
ing to decide which solutions of the EKG system are of 
physical interest, it also simplifies the numerical solu- 
tions. 

The exponential case fits well into the dynamical sys- 
tem formalism, as has been extensively shown in the liter- 
ature, but the formalism can be easily adapted for other 
scalar potentials. We expect to report on them elsewhere. 
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